A class of quantum many-body states that can be efficiently simulated 
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We introduce the multi-scale entanglement renormalization ansatz (MERA), an efficient represen- 
tation of certain quantum many-body states on a D-dimensional lattice. Equivalent to a quantum 
circuit with logarithmic depth and distinctive causal structure, the MERA allows for an exact evalu- 
ation of local expectation values. It is also the structure underlying entanglement renormalization, a 
coarse-graining scheme for quantum systems on a lattice that is focused on preserving entanglement. 



PACS numbers: 

A better understanding of quantum entanglement has 
enabled significant progress in the numerical simulation 
of quantum many-body systems over the last few years 0, 
E3- 01 ■ Buildingon the density matrix renormalization 
group (DMRG) 0-a well-established technique for one- 
dimensional (ID) systems on a lattice-, new insight from 
quantum information science has led, e.g., to efficient 
algorithms to simulate time-evolution [l| and address 2D 
systems 01 • 

A key ingredient of such algorithms is the use of a net- 
work of tensors to efficiently represent quantum many- 
body states. Examples of tensor networks include matrix 
product states (MPS) for ID systems, tree tensor net- 
works (TTN) jj] for systems with tree shape, and pro- 
jected entangled-pair states (PEPS) Q for 2D systems 
and beyond, the three structures differing in the graph 
that defines how the tensors are interconnected into a 
network: the graphs for MPS, TTN and 2D PEPS are, re- 
spectively, a chain, a tree and a 2D lattice. Importantly, 
from these tensor networks the expectation value of lo- 
cal observables can be computed efficiently. But whereas 
from a MPS and a TTN such calculations are exact, from 
a PEPS -which has a much wider range of applications- 
local expectation values can only be obtained efficiently 
after a number of approximations. 

In this manuscript we present the multi-scale entan- 
glement renormalization ansatz (MERA), a structure 
that efficiently encodes quantum many-body states of 
D-dimensional lattice systems and from which local ex- 
pectation values can be computed exactly. A MERA 
consists of a network of isometric tensors in D+l dimen- 
sions, where the extra dimension can be interpreted in 
two alternative ways: either as the time of a peculiar 
class of quantum computations, or as parameterizing dif- 
ferent length scales in the system, according to succes- 
sive applications of a lattice coarse-graining procedure 
known as entanglement renormalization Q. Focused on 
preserving entanglement, MERA is a promising candi- 
date to describe emergent quantum phenomena, includ- 
ing quantum phase transitions, quasi-particle excitations 
and topological order. Here we establish its connection 
to entanglement renormalization and explore some of its 
basic properties: efficient contractibility of the network, 




FIG. 1: Quantum circuit C transforms state |0) into the 
A-site state |^) of a ID lattice C. C contains 2 A — 1 gates 
organized in 0(log A) layers labeled by a discrete time 9. 



leading to an efficient evaluation of local expectation val- 
ues; inherent support of algebraically decaying correla- 
tions and of an area law for entanglement; versatility to 
adapt to both the local and global structure of the sys- 
tem's underlying lattice; and ability to naturally assimi- 
late symmetries such as invariance under translations or 
rescaling - which result in substantial gains in computa- 
tional efficiency. 

Let us consider a square lattice C in D spatial dimen- 
sions consisting of N sites, where each site s E C is de- 
scribed by a complex vector space V of finite dimension 
X- Let I'i') S V® N denote a pure state of lattice C and 
= tr^ ( | \D r ){\E f | ) its reduced density matrix for site s. 
We study states |\&) that can be generated by means of a 
certain quantum circuit C of depth 9 = 2 log 2 (N) — 1, see 
Fig. for the D — 1 case. Each site in C corresponds 
to one outgoing wire of C, and we label both with the 
same index s. 

The causal cone C^l of outgoing wire s, defined as the 
set of gates and wires that may influence state pM, plays 
an important role in our discussion. The key feature of 
circuit C is that each causal cone has bounded width, 
that is, the number of wires in a time slice C^ s ' of is 
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FIG. 2: A MERA M inherits the causal structure of quan- 
tum circuit C: the causal cone C' s ' for site s has bounded 
width. When reversing the arrow of time 9, M implements 
entanglement renormalization transformations. 

bounded by a constant that is independent of N. [Here 
9 labels a discrete time, with 9 £ {0, l,--- ,6}]. For 
instance, the causal cone of any outgoing wire in Fig. (JTJ 
involves at most four wires at any time 9, as highlighted 
in Fig. @. More generally, the causal cone of outgoing 
wires in the D-dimensional case can be seen to comprise 
at most 3 D_1 4 wires at a given time 9, see Fig. (JSJl. 

A MERA for |W) is a tensor network M. that corre- 
sponds to quantum circuit C except for some minor cos- 
metic changes. Each two-body unitary gate u of C gives 
rise to a tensor of M.. But incoming wires in state |0) are 
eliminated in the tensor network, producing three kinds 
of tensors: (i) the top tensor t — u\Q}\0) of M. has two 
indices and is normalized to 1, 

(*W = («)#k/s=o, E(*V(V = 1 ; (!) 

(JLV 

(ii) tensors in every second row are isometries w = u\0), 

= «£i/3= , = <w; ( 2 ) 

(%%%) the rest of tensors in A4 are unitary gates ti, 

Y^W'gM'tf =&* a >6w, (3) 

EK)^w;t=v^, (4) 

a/3 

that we call disentanglers for reasons that will become 
clear later. Notice that the computational space required 
to store M grows as 0(x A N), that is, linearly in N, given 
that there are 2N — 1 tensors and each tensor depends 
on at most % 4 parameters. 

Thus a MERA is an efficient representation of \^) con- 
sisting of a tensor network M in D + 1 dimensions with 
two properties: (i) tensors are constrained by Eqs. 




FIG. 3: Detail of a causal cone in a 2D MERA, implying that 
its width is bounded. 3x3 sites (bottom left) are mapped 
into 3x3 sites (top right) by layers of tensors that act along 
the y and x direction [some stages involve 3x4 sites]. 

(J2J or J2J-I0J; (ii) each open wire s, associated to one site 
of the underlying lattice £, has a causal cone with 
bounded width. As a consequence of this peculiar causal 
structure, the reduced density matrix of a small number 
of lattice sites can be computed exactly with remarkably 
small cost. In what follows p\ and P2 are integers that 
depend on the spatial dimension D of C. 

Lemma 1: The one-site density matrix pW can be 
computed from a MERA with time 0(x Pl log AT). 

Proof: For each time slice of the causal cone C^ s \ 
we compute its reduced density matrix o~g. As shown in 
fig. (0J for D = 1, ag+i can be obtained from og with 
a cost polynomial in \ an d independent of N . Recall, 
finally, that pM = cr and that 9 = 2\ag 2 (N) - 1. 

Lemma 2: The two-site density matrix p^- SlS2 ^ can be 
computed with 0(x P2 log AT) time. 

Proof: Again, the causal cone Cl slS2 ] has logarithmic 
depth and a width independent of N, see Fig. (0, and 
the reduced density matrix o~g+\ for time slice C^_t* can 

be obtained from ag for df 1 ^ at bounded cost. 

The above lemmas can be easily extended to a /c-site re- 
duced density matrix p[ Sl *'" s w (with computational time 
scaling exponentially in k) and they imply that the ex- 
pectation value of local observables, such as two-site cor- 
relators 

C 2 {s 1 ,s 2 ) = (^\A^B^\t>) =tiip [siS2] A^B^), 

(5) 

can be computed efficiently. 

We have defined a MERA for |W) in terms of a quan- 
tum circuit C that transforms a product state |0) 8Ar into 
|$) by means of layers of unitary gates. Let |^ T ) 
denote [the non-trivial part of] the state of circuit C at 
time 9 = — 2r, with j^o) = \^)- An alternative inter- 
pretation of the MERA can be obtained by considering 
the sequence of states {I'Fo), 1^2), ■ ■ ■ }, which cor- 
respond to undoing the quantum evolution of C back in 
time. Notice that |\& r +i) is obtained from \^> T ) by ap- 
plying two layers of tensors in A4 . The first layer is made 
of disentanglers that transform |'F T ) into a less entangled 
state \^' T ). The second layer is made of isometries that 
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FIG. 4: Efficient computation of p^. (i) Part of the causal 
cone C' s ' of p' s ', D = 1. (ii) How to compute o"2fc+i from 
<T2fc_i through simple tensor multiplications (Hi) and (it)). 

combine pairs of nearest neighbor wires into single wires, 
turning the state \^>' T ) of N T wires into the state l^r+i) of 
N T /2 wires, where N T = 2 e ~ r . That is, M. implements a 
class of real space coarse-graining transformations known 
as entanglement renormalization |7j. 

More generally, we can use the MERA to transform the 
sites of lattice Co = C, as well as operators defined on C 
such as a Hamiltonian Hq, so as to obtain a sequence of 
increasingly coarse-grained lattices {Co, C\, £2, ■ ■ • } and 
corresponding effective Hamiltonians {Ho, H±, if 2, ■ ■ ■ }• 
Since an operator defined on site s G C is mapped 
into an operator contained in the causal cone CW, lo- 
cal operators in C remain local when mapped into C T , 
1 < t < log 2 (A) — 1. Notice that one site in C T is ob- 
tained from a D-dimensional hypercube with 0(2 T ) sites 
in C. Thus |\I/ T ) can be understood as retaining the struc- 
ture of at length scales greater than T~I D . 

In it has been shown that a MERA can encode, in 
a markedly more efficient way than a MPS ||, accurate 
approximations to the ground state of a quantum criti- 
cal system in a ID lattice. There correlators C<x{s\,S2) 
decay as a power law with the distance r between sites 
s\ and S2 and the entanglement entropy S(p^ B '^) of a 
block Bi of I adjacent sites scales as logZ 0- Next we 
give an intuitive justification why MERA supports alge- 
braic decay of correlations in any dimension D. We also 
show that a MERA can carry an amount of block entan- 
glement compatible with logarithmic scaling in ID and 
with a boundary law Z D_1 for D > 1. 

It follows from the A-shaped causal cone & lS ^ that 
we can compute p[ slS2 l from the density matrix p f cor- 
responding to a hypercube of (at most) 4 D sites of Cf, 
where f sa D log 2 r, by means of a sequence of density 
matrices {p f , . . . , p 2 , Pi, p' SlS2 ' }, see Fig. That is, 

p[sis 2 ] j g obtained from p> after O(logr) transformations. 




FIG. 5: Causal cones C [siS2] and C [Bl] for a two-site density 
matrix p' slS2 l (left) and an /-site density matrix p' Sl ' (right). 

If, as it can be argued in a (scale invariant) critical ground 
state, each of these transformations reduces correlations 
by a constant factor z < 1 , we readily obtain a power law 
scaling for 6*2(51, S2), 

C 2 (si, S 2)-z losr = r^, g = logi. (6) 

This property is enabled by the fact that two sites at a 
distance r in C are connected through a path of length 
O(logr) in A4, and it indicates that a MERA is par- 
ticularly suited to describe states with quasi-long-range 
order, such as critical ground states. 

We now consider the entropy S(p^- Bl ^) for the density 
matrix p^ B '^ of a hypercube Bi made of l D sites. The 
causal cone of J5/ shrinks exponentially fast with r and 
we see that, once more, we can compute p>- Bl ' from the 
density matrix rj^ for a hypercube made of (at most) 
4 sites of lattice Cf, for f w D\og 2 l, through a se- 
quence of density matrices {r]f , • • ■ ,r]i,p^ B '^}, see Fig. 
(|5|). Density matrix rj T is obtained from rj T+ i by: (i) 
applying a layer of isometries and a layer of disentan- 
glers that do not change the entropy of ?7 T +i; (ii) trac- 
ing out n T — 0{2^ T ~ T ^ D ^ 1 ^ D ) boundary sites. Tracing 
out one boundary site increases the entropy by at most 
log 2 x bits, the total increase in entropy AS T being at 
most n T log 2 (x) bits. Thus the entropy of p^ B '^ fulfils 

f r 

S(p^) - S( Vf ) = < log 2 ( X )5> T , (7) 

T=l T=l 

where S(r]f) is at most 4 D log 2 x- 

In a ID MERA the number n T of sites that are traced 
out is bounded by a constant c [i.e., a hypercube has only 
two boundary sites] for any r, 1 < r < f = O(log^), and 

S( P{1} ) - S( Vf ) < log a (x)cf - OQoel). (8) 

As numerically confirmed in [viEr,evER], in a MERA for 
critical ID ground states AS T is independent of r, leading 
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to the logarithmic scaling S(p^- B '^) w log(Z), whereas for 
a non-critical ground state A5 T vanishes for r 3> log 2 £, 
where £ is the correlation length, so that S(p^ B '^) satu- 
rates for I 3> £. 

In Z) > 1, instead, the n T decays exponentially with r, 
the upper bound for the «S(pI B, J) being dominated by the 
contribution from small t, 

T 

S(p^)-5(Tfr) < log 2 ( x )^2^)^ 

T=l 

« log 2 ( X ) 2^ » log 2 ( X ) 2^ lo ^W^ 

= log 2 ( x ) Z^" 1 . (9) 

That is, a MERA in D > 1 supports block entanglement 
that scales at most according to a boundary law S^- B '^ « 

For the sake of concreteness, we have analyzed the 
case where £ is a square lattice. The structure of the 
MERA, however, can be adapted to a more generic lat- 
tice, with arbitrary local, geometric and topological prop- 
erties, while preserving its distinctive causal structure. 
For instance, in D = 2 dimensions a specific MERA can 
be built to represent states of a triangular lattice; or of 
a lattice with random vacancies or linear or bulk defects; 
or to account for a variety of boundary conditions (e.g., 
plane, cylinder, sphere or torus). In addition, the num- 
ber of levels x can vary throughout M.. Adjusting the 
MERA to the specifics of a problem often leads to com- 
putational gains. 

In particular, the symmetries of |^} can be assimilated 
into the MERA. An internal symmetry, such as SU(2) in- 
variance, results in a series of constraints for the tensor 
in M which then depend on less parameters. More 
substantial gains are obtained when is invariant un- 
der translations [that is, cyclic shifts by one lattice site 
in a system with periodic boundary conditions] , since all 
the tensors in a layer of M. can be chosen to be the same 
and the MERA depends only on about \ 4 l°g2 N param- 
eters. But the most dramatic savings occur for states 
that are invariant under entanglement renormalization 
transformations, even in an infinite lattice Here 
all tensors in M. are the same and the MERA depends 
on just Of x ) parameters. Scale invariant critical ground 
states [lof can be shown to belong to this class. 

We conclude with a few pointers to future work. On 
the one hand, most techniques to simulate quantum sys- 
tems with a MPS can be generalized to a MERA. This 
include algorithms to compute the ground and thermal 
states, and to simulate time evolution. For instance, in 
order to update M. after the action of a unitary gate 
y[ SlS2 ] on sites s± and S2, only the tensors in the causal 
cone Cl SlS2 l need to be modified, with a cost logarithmic 
in N. 

On the other hand, the potential of a MERA is not re- 
stricted to the representation of individual states. Notice 



that by feeding the incoming wires, labeled by r, of quan- 
tum circuit C with state <8>£Li |<4- r ') instead of ®£L]J0M), 
we can generate an infinite family of entangled states 
au represented by a MERA that only differs in 
the isometric tensors and such that, given the unitarity 
of C, fulfill 

<*{Mi*w}Hn^w>. (io) 

r 

This can be used to encode, in just one (generalized) 
MERA, not only the ground state of a lattice system but 
also its quasi-particle excitations. This is also used, in 
systems with topological quantum order, to store topo- 
logical information in the top tensor of a MERA ^1 ■ 
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